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ABSTRACT 



Context. We present the AGILE gamma-ray observations of the field containing the puzzling gamma-ray source 3EG 11835+5918. 
This source is one of the most remarkable unidentified EGRET sources. 

Aims. An unprecedentedly long AGILE monitoring of this source yields important information on the positional error box, flux 
evolution, and spectrum. 

Methods. 3EG J1835+5918 has been in the AGILE field of view several times in 2007 and 2008 for a total observing time of 138 
days from 2007 Sept 04 to 2008 June 30 encompassing several weeks of continuous coverage. 

Results. With an exposure time approximately twice that of EGRET, AGILE confirms the existence of a prominent gamma-ray source 
(AGL J 1836+5926) at a position consistent with that of EGRET, although with a remarkably lower average flux value for photon 
energies greater than 100 MeV A 5-day bin temporal analysis of the whole data set of AGL J1836+5926 shows some evidence for 
variability of the gamma-ray flux. The source spectrum between 100 MeV and 1 GeV can be fitted with a power law with photon 
index in the range 1.6-1.7, fully consistent with the EGRET value. 

Conclusions. The faint X-ray source RX J1836.2+5925 that has been proposed as a possible counterpart of 3EG J1835+5918 is well 
within the AGILE error box. Future continuous monitoring (both by AGILE and GLAST) is needed to confirm the gamma-ray flux 
variability and to unveil the source origin, a subject that is currently being pursued through a multiwavelength search for counterparts. 

Key words, gamma-rays: observation; 



1. Introduction 

The gamma-ray source 3EG J1835+5918 (also known as GEV 
J1835+5921 and GRO J1837+59) is one of the most puz- 
zling high-energy sources in the sky, and it has been the sub- 
j ect of considerab le interest since its discovery by EGRET 
dNolan et al. 19941) . This source was catalogued among the 
brightest high-latitude unidentified sources with an average 
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gamma-ray flux above 100 MeV of ^ egret - (60 + 4) 10"** 
photons cm"^ s" ' , and a hard gamma- ray spectrum of photon in- 
dex y - 1.69 + 0.07 (iHartman et al. 1 999). Owing to the spec- 
tral hardness (similar to that of the Vela pulsar), the source flux 
near 1 GeV is expected to be about half of the Crab value. The 
EGRET error box is well-determined, resulting in an error radius 
of 8'at the 95% confidence level. 

The source 3EG J1835+5918 is located 25° off' the Galactic 
plane and therefore is not significantly affected by the Galactic 
diffuse gamma-ray emission. EGRET pointed at this high- 
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galactic latitude region several times ( iHartman et al. 19991) . 
The resulting position, gamma-ray lightcurve, variabil- 
ity, and spectrum have been discussed in several papers : 
dMcLaughlin et al 1996l:lNolan et al. 19"96HReimer et al. 20011) . 
The analysis of the whole EGRET dataset produces a 
gamma-ray flux database consistent with being constant 
(IHartman et al. 19991: iReimer et aL 200lh . However, it is impor- 
tant to notice that a claim for variability of 3EG J1835+5918 
was indicated by Nolan et al. (1996). 

The relatively strong gamma-ray flux, the hard spectrum, and 
the lack of an obvious blazar counterpart in the EGRET error 
box of 3EG J1835+5918 spurred considerable interest from 
several observing groups. After extensive X-rays, rad io, and 
optical coverage o f the EGRET error box (Mirabal et a L2000l: 
iReimer et aL 20011: iTotani et al. 20021: iHalpern et al. 20021) the 
search for possible counterparts singled out RX J1836.2-I-5925, 
a relatively faint soft X-ray source with no radio or optical emis- 
sion: a set of characteri stics similar to thos e of middle-aged 
radio-quiet neutron stars dCaraveo et al. 1996h . Thus, on the ba- 
sis of the gamma-ray properties of 3EG J1835+5918, coupled 
with those of RX J1836.2-H5925, the sources was named the 
next Geminga by ( Halpern et al. 2002|) owing t o its Geminga- 
like phenomenology ( Bignami and Caraveo 1996) . A search for 
X-ray pulsation, with period values ranging from 1 ms to 10 
s, yielded a 35% u pper limit to the source-pulsed fraction 
jHalpern et al. 2007h . In view of the much lower pulsed fractions 
measured for all known middle-aged neutron stars, such a result 
did not weaken this tentative identification. 

Taking advantage of the AGILE unprecedented exposure 
gathered during the science verification phase and early point- 
ings of the Cycle- 1 program, in this Letter we are now able to 
assess the source characteristics anew. Section 2 describes the 
AGILE-GRID observations and main results, while their impli- 
cations are discussed in Section 3. 

2. The AGILE GRID observations of 
3EGJ1 835+591 8 

The AGILE istrument jTavani et al. 20081) is composed 
of three detectors: a Tungsten-Si licon Tracker (ST) 
jBarbielfini et al. 20021: iPrest et al. 2003h . with a large field 
of view (~ 2.5 sr), optimal time resolution and angular reso- 
lution, and good sensitivit y; a Silicon- based X-ray detector, 
Super-AGILE (SA) (Feroci et al. 2007h for imaging in the 
energy ran ge 18 keV - 60 keV and a CsI(Tl) Mini-Calorimeter 
(MCAL) dLabanti et al. 20061) that detects gamma-rays or 
particle energy depositions between 300 keV and 100 MeV. 
The ST and MCAL form the AGILE Gamma-Ray Imaging 
Detector (GRID) for observations in the gamma-ray energy 
range 30 MeV - 50 GeV. The instrument is com pleted by an 
anti-coincidence (AC) system (j Perotti et al. 20 0^. made with 
plastic scintillator layers, for the rejection of charged particles, 
and by an e fficient trigger log ic for gamma-ray and X-ray data 
acquisition jTavani et al. 2008 !). 

Data analysis was performed with the BUILD 15 of the 
AGILE Standard Pipeline, publicly available at the ASI Data 
Center web site (http://agile.asdc.asi.it/). The data were pro- 
duced starting from the Level- 1 data. The events collected dur- 
ing the passage in the South-Atlantic Anomaly and the Earth 
albedo background were consistently rejected. The GRID event 
direction were reconstructed by a Kalman filter technique. To 
reduce the particle background contamination, we selected with 
the GRID filter FT3ab^ only events flagged as confirmed y- 
ray events (G class events, corresponding to a sensitive area of 



~ 300 cm^ at 100 MeV). All the fluxes reported in this Letter 
have been consistently checked and scaled with their corre- 
sponding errors using the ratio R - '^'!^y'ja)iaLE '^o^^P'^ted for 
similar off-axis angles, with <^{Veld) the gamma-ray flux above 
100 MeV reported by EGRET and measured by AGILE. Scaling 
our fluxes measured for 3EG J1835+5918 to the Vela PSR mea- 
surements by AGILE and EGRET also takes systematic effects 
into account (estimated to be of the order of the 10% of the fluxes 
reported here). 

All AGILE observations of the field containing the 
3EG J1835H-5918 are listed in Table [T] The columns report the 
viewing periods name (as reported in the ASI Data Center web 
site), the galactic coordinates of the centre of the field of view, 
the start and end of the observation in UTC, the average source 
off-axis angle, the period-averaged flux Fy (E > 100 MeV) in 
photons cm"^ s '(if y/TS < 3 a 2 - cr upper limit is reported), 
the I - cr statistical and systematic uncertainty of the flux, the 
number of counts, and in the last column the statistical signifi- 
cance of the likelihood ratio test. 

AGILE counts, exposure, and galactic background maps 
were generated with a bin size of 0.25° x 0.25° for E > 
100 MeV to compute the source flux and its evolution, 
while maps with a bin size of 0.1° x 0.1° were used to op- 
timi ze the source p osition using a likelihood analysis method 
(Ma ttox et al. 1996a) . Our analysis was performed over a region 
of 10° radius. 

2.1. The gamma-ray source positioning 

To determine the most likely location of the gamma-ray source, 
we used only the best viewing periods, i.e. those in which the 
source was within 50° from the instrument pointing direction. 
For E > 100 MeV, the best position computed by the like- 
lihood analysis is / = 88.90° b - 24.99°, corresponding to 
a = 18°36'20".ll, 6 = 59°26'42".3. We name this source 
AGL J1836-I-5926. Figure [T] shows a Gaussian-smoothed counts 
map of the source with our 95% statistical and systematic error 
ellipse with semiaxis a = 0.265° and b = 0.200°. To conser- 
vatively take both systematic and statistical effects into account, 
we obtain a 95% confidence level radius by linearly adding a 
systematic error of 0.1° to the 95% statistical error ellipse. 

The AGILE 95% e rror ellipse is consis tent with the error box 
of 3EG J1835-H5918 (Matt ox et al. 20011) . and the distance be- 
tween AGILE and EGRET centroids is ~ 10.2 '. The distance be- 
tween the position of the X-ray source RX J1836. 2-1-5925 and the 
AGILE 95% maximum likelihood contour level barycentre is ~ 
1 .4 '. The source detection significance, as derived from a maxi- 
mum likelihood analysis, is 16.7-cr. A total of 499 photons with 
£■ > 100 MeV have been collected. The source flux, averaged 
over the whole AGILE dataset, is (^acile = (38.7 + 3.0) x 10"** 
photons cm"^ s"' , which turns out to be approximately 60% less 
than the value reported in the 3rd EGRET Catalogue. 

2.2. The gamma-ray spectrum 

Our standard GRID spectral analysis was performed using 3 en- 
ergy bins, namely 100-200 MeV, 200-400 MeV, 400-1000 MeV. 
Photons below 100 MeV and above 1 GeV have not been used 
in the analysis reported here. The analysis was obtained by sum- 
ming the observation periods OB5600, OB5700, and OB5800 
(see Table[Tli. A power-law model fit to the data yields a photon 
index y = 1.65 ±0.22. 
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Table 1. AGILE observations ofSEG J1835+5918. 



Viewing period 


1 


b 


start (U 1 C) 


End (U IC) 


Off-axis angle 




AF-y 


Counts 




OB2300 


134.882 


11.822 


'~\c\f\^ c\f\ r\ A 1 f\f\ 

2007-09-04 12:00 


2007-09-12 12:00 


50 


< 40 




< 13 


1.8 


OB4610 + OB4630 


143.365 


41.588 


^f\f\^ f\ ^ A f\Ci f\f\ 

2007-10-24 08:00 


2007-11-01 12:00 


50 


81 


18 


37 ± 10 


4.4 


OB4800 


69.594 


4.623 


2007-11-02 12:00 


2007-12-01 12:00 


27 


40 


6 


99 ± 16 


7.7 




OO.Ol J 




ZUU/-1Z-U1 iz.uu 


9007 1 7 O^ OQ-OO 

zuu/-iz-uj uy.uu 




jy 




18-1-^ 
lo ± J 


.:>.o 


OB4910 + OB4920 


85.119 


-9.416 


2007-12-05 09:00 


2007-12-16 12:00 


38 


48 


12 


44 ± 10 


5.2 


OB5210 


77.309 


40.628 


2008-02-09 09:00 


2008-02-12 12:00 


20 


60 


19 


21 ± 7 


4.2 


OB5600 


53.039 


6.474 


2008-04-10 12:00 


2008-04-30 12:00 


40 


40 


10 


57 ± 13 


5.0 


OB5700 


104.852 


35.439 


2008-04-30 12:00 


2008-05-10 12:00 


30 


49 


11 


49 ± 10 


6.2 


OB 5 800 


74.05 


0.273 


2008-05-10 12:00 


2008-06-09 12:00 


30 


50 


9 


129 ± 18 


9.6 


OB5820 


93.6 


-1.16 


2008-06-15 12:00 


2008-06-30 12:00 


30 


25 


8 


37 ± 11 


4.2 




Fig. 1. Gaussian-smoothed counts map with squared scale of 
AGL J1836+5926 with photons of E > 100 MeV. The red 
curve is the EGRET 95% error ellipse, the blue curve is the 
AGILE 95% statistical plus systematic error ellipse, the cross 
is the AGILE 95% error ellipse barycentre, and the black circle 
is the position of the X-ray source RX J1836. 2+5925. 



2.3. The gamma-ray lightcurve and variability analysis 

The search for the variabiHty of the gamma-ray emission of 
3EG J1835H-5918 produced somewhat inconclusive results. 
Using the first EGRET viewing periods, McLaughlin et al. 
(1996) and Nolan et al. (1996) classified 3EG J1835+5918 as 
a possibly variable source. On the other hand, using the whole 
EGRET data set Reimer et al. (2001) and Nolan et al. (2003) did 
not confirm the previous hint of variability. (We note that none 
of the previous EGRET variability studies used the same 5-day 
intervals adopted in this Letter) 

To study the source temporal behaviour, each AGILE view- 
ing period was first analysed to determine the parameters of 
the Galactic diffuse radiation model and of the isotropic dif- 
fuse intensity level. The source flux density was then estimated 
independently for each 5-day temporal bin with the gamma- 
ray model parameters fixed at the values obtained in the first 
step. Figure |2] shows the AGL J1836-t-5926 gamma-ray flux 
lightcurve during the period September 2007 to June 2008 for 
photons with £ > 100 MeV. 



We fit the AGILE GRID fluxes to a constant value 
(the weighted mean of the average flux values), and com- 
puted the_J^_j[ariabilitY_coefficient for the AGILE observa- 
tions jMcLaughlin et al. 1996) . The 2-cr upper limits (UL) were 
properly treated, a ssigning a value and a cr equal to UL/2 
(Tor res et al. 2001a '). A value of V > 1 is an indication of vari- 
ability within the observing period. The AGILE observation 
block durations of Table [1] span a non-homogeneous range of 
timescales from 1 week to more than a month. To study the flux 
variability of AGL J1836+5926 on a timescale of a few days, 
we chose a temporal bin size of 5 days as a good compromise 
between statistics and variability characteristics. Using a 5-day 
bin size, we 

obtain ^ 38.34 for 24 de grees of freedom and a 
value of y = 1.495. We can exclude the flux being constant at 
the 96.80% level. For a temporal bin of 4 days we obtain V ^ 1. 
For temporal bins of 3 days or shorter, there are too many up- 
per limits during the whole exposure time. For temporal bins of 
more than 5 days, we obtain V < 1 , as expected for averaging a 
signal potentially variable on a timescale of a few days. 

We note that, during the long uninterrupted coverage ob- 
tained in Nov-Dec 2007 and in April- June 2008, the gamma-ray 
source AGL J 1836+5926 was not detectable (failing to reach 
the 3-cr threshold) for about half of the observing time. We also 
determined hard X-ray upper limits for AGL J1836-t-5926 in 
the energy range 18-60 keV from an analysis of Super- AGILE 
data. A 3-sigma upper limit of 10 mCrab was obtained for the 
OB5700. 



2.4. EGRET vs AGILE 

After one year of operations, AGILE has doubled the EGRET 
exposure time of 3EG J1835-I-5918, confirming the source po- 
sition and spectral shape and showing evidence of variabil- 
ity. A comparison of the AGILE and EGRET observations 
of 3EG J1835H-5918 is reported in Table ^ With a few ex- 
ceptions, the EGRE T pointings were generally one-week long 
(iReimer et al. 200 ih . and they were not optimally suited to trac- 
ing an emission variable on week to month long timescales. On 
the other hand, AGL J1836-H5926 has been continuously in the 
AGILE field of view for more than two months in April-June 
2008, in addition to several weeks scattered between September 
2007 and December 2007. It is too early to say whether the 
apparent variations are erratic or follow some kind of cyclic 
pattern. Longer, continuous exposures (both with AGILE and 
GLAST) are needed to clarify this crucial point. 
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Fig. 2. AGILE gamma-my light-curve ofAGL J1836+5926 for photons above 100 MeV with a 5-day temporal bin. 



Table 2. Comparison of total exposure and number of photons 
gathered by EGRET and AGILE from 3EG J 1835+591 8 for E > 
100 MeV. 





cm^ s sr/Q.5° bin 


Photons £ > 100 MeV 


EGRET 


52900 


452 


AGILE 


118325 


499 



3. Discussion and conclusions 

The AGILE observations and monitoring of AGL J 1836-1-5926 
adds relevant information about this puzzling source. Postponing 
an account of our multifrequency observations of the region to 
a forthcoming paper, we briefly outline here a few important 
points regarding the search for a counterpart. We confirm a point 
already noticed by other authors, i.e., the absence of a blazar or 
any other relatively bright radio sources in the field containing 
the AGILE (and EGI^ET) error box. The analysis of Mattox et al. 
(2001) shows a 100 mJy source (B1834-H5904) positioned at the 
99.5% probability contour (12'. 7 from the EGRET centroid po- 
sition and 2 1 ' .8 from the AGILE centroid position). By using the 
NED, SIMBAD, and the Massaro et al. (2007) catalogues, we 
confirm the absence of a radio-loud blazar in the revised AGILE 
error box. 

Outside this error box, we notice the existence of RGB 
J1841H-591, a BL Lac type object positioned at 43'.5 from the 
AGILE centroid position. Even though occasional contributions 
from this object cannot be excluded in our analysis of AGILE , 
we emphasise that the AGILE integrated flux positioning at the 
95% contour level is clearly not consistent with any substantial 
contribution from RGB J 1 84 1 -1-59 1 . 

We also notice fliat the faint X-ray source RX J1836.2-H5925 
that was proposed as a possible counterpart is well within the 
AGILE error box of AGL J 1836+5926. A hint of variability of 
this source was noticed when comparing t wo HRI ROSAT ob - 
servations taken almost three years apart dMirabal et al. 200T!) . 



prompting a claim f or the discovery of a new class of compact 
gamma-ray sources dMirabal et al. 2000h . However, subsequent 
Chandra X-ray observations showed a practically constant X-ray 
flux of RX J 1836.2+5925, and weakened the variability claim 
in favour of a scenario encompassing a constant Geminga-hke 
source. 

We note here that our group observed the region contain- 
ing RX J1836.2+5925 with the SWIFT XRT and XMM-Newton 
during the period May -June 2008. Analysis of these data and a 
full discussion of the AGL J 1836+5926 counterpart problem 
will be presented elsewhere. 

Future continuous monitoring (both by AGILE and GLAST) 
is needed to confirm the gamma-ray flux variability and to unveil 
the source origin. 
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